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Abstract — We establish a lower bound on the noncoherent 
capacity pre-log of a temporally correlated Rayleigh block-fading 
single-input multiple-output (SIMO) channel. Our result holds 
for arbitrary rank Q of the channel correlation matrix, arbitrary 
block-length L > Q, and arbitrary number of receive antennas R, 
and includes the result in Morgenshtern et al. (2010) as a special 
case. It is well known that the capacity pre-log for this channel 
in the single-input single-output (SISO) case is given by 1 — Q/L, 
where Q/L is the penalty incurred by channel uncertainty. Our 
result reveals that this penalty can be reduced to 1 /L by adding 
only one receive antenna, provided that L > 2Q 1 and the 
channel correlation matrix satisfies mild technical conditions. 
The main technical tool used to prove our result is Hironaka's 
celebrated theorem on resolution of singularities in algebraic 
geometry. 

I. Introduction 

It was shown in HI that the noncoherent capacitjQ pre- 
log for single-input multiple-output (SIMO) correlated block- 
fading channels can be larger than the pre-log in the single- 
input single-output (SISO) case. This result was surprising as 
it disproved a conjecture in ||2l on the pre-log in the SIMO 
case being the same as in the SISO case. The channel model 
analyzed in fTl assumes that the fading process is independent 
across blocks of length L and temporally correlated within 
blocks, with the rank of the corresponding L x L channel 
correlation matrix given by Q < L and the number of receive 
antennas R = Q. For this channel model, under additional 
technical conditions on the channel correlation matrix, a pre- 
log of (1 — 1/L) was established in fll. In contrast, in the 
SISO case the pre-log is given by \ — Q/L. 

The assumption R = Q made in |[T] is very restrictive, 
and the proof technique used in |T] heavily relies on this 
assumption. More precisely, the main result in H] is based 
on a lower bound on the differential entropy h{y) of the 
channel output signal that is obtained by applying a change of 
variables argument [[T] Lem. 3]. The proof is then completed 

The work of Erwin Riegler was partially supported by the WWTF project 
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' Noncoherent capacity denotes capacity in the setting where transmitter 
and receiver know the channel statistics but neither of them is aware of the 
channel realizations. 



by showing that the expected logarithm of the Jacobian de- 
terminant corresponding to this change of variables is finite. 
For R < Q, the Jacobian determinant takes a very involved 
form making it difficult to say anything about its expected 
logarithm. The main contribution of this paper is to resolve 
this problem by introducing a new proof technique based on a 
result from algebraic geometry, namely ||3] Th. 2.3], which is 
a consequence of Hironaka's celebrated theorem on resolution 
of singularities JU, ||5]- Roughly speaking, this result allows 
us to rewrite any real analytic function ||6] Def. 1.6.1] locally 
as a product of a monomial and a nonvanishing real analytic 
function. The proof of our main result, a lower bound on the 
pre-log for the correlated block-fading channel with arbitrary 
number of receive antennas R < Q, is then effected by using 
this factorization to show that the integral of the logarithm of 
the absolute value of a real analytic function over a compact 
set is finite, provided that the real analytic function is not 
identically zero. This method is very general and could be 
of independent interest when one tries to show that a certain 
differential entropy is finite. 

We conclude by noting that the main result in this paper 
shows that the pre-log penalty Q /L incurred in the SISO case, 
which is due to channel uncertainty, can be reduced to 1/L 
by adding only one receive antenna (i.e., by taking R ~ 2), 
provided that L > 2Q — 1 and the channel correlation matrix 
satisfies mild technical conditions. In the limit L,Q ^ oo with 
L/Q constant, the 1/L penalty in the SIMO case becomes 
arbitrarily small, whereas the Q/L penalty in the SISO case 
remains unchanged. 

A. Notation 

Finite subsets of the set of natural numbers, I C N, 
are denoted by calligraphic letters and we write \I\ for the 
cardinality of I. We use [m : n] to designate the set of natural 
numbers {m, m+1, . . . , n}. Uppercase boldface letters denote 
matrices, lowercase boldface letters designate vectors. The 
superscripts ^ and ^ stand for transposition and Hermitian 
transposition, respectively. The all-zero matrix of appropriate 
size is written as 0. For a matrix A e C*^^^, the entry in the 



ith row and jth column is denoted by ai,j and we write aj 
for its ith row. If I C [1 : Af], we denote by Ai G C'^l the 
submatrix of A obtained by retaining all rows aJ of A with 
row index i E I. Similarly, for a vector x E C^^ we denote 
its ith entry by a;^ If I C [1 : M] we denote by Xj G C'-^' the 
vector obtained by retaining the entries Xi of x with i G I. We 
write Ci for the ith unit vector of appropriate size and Im for 
the identity matrix of size AI x AI. For a vector x, diag(a;) 
denotes the diagonal matrix that has the entries of x in its 
main diagonal. For two matrices A and B of arbitrary size, 
diag(A, B) is the 2x2 block matrix that has the matrix A as 
upper left block, B as lower right block, and as upper right 
and lower left block. For N matrices Ai, . . . , Ajv, we define 
diag{Ai,. . . ,An) = diag(diag(Ai, . . . , Aat^i), Aat). We 
designate the Kronecker product of the matrices A and B 
as A (E) B; to simplify notation, we use the convention that 
the ordinary matrix product precedes the Kronecker product, 
i.e., AB C = (AB) (g) C. For a function /, we write 
/ ^ if there exists a vector x in the domain of / such that 
f{x) 7^ 0. For two functions / and g, the notation / = 0{g) 
means that limsup^^^(|/(x-)/5(a;)|) < od. If x G K, [x\ ^ 
max{?n E Z \ m < x} and [x] = min{m £ X \ m > x} 
with Z denoting the set of integers. The logarithm to the base 
2 is written as log(-). The expectation operator is denoted by 
E[-]. Finally, CAf{a, C) stands for the distribution of a jointly 
proper Gaussian random vector with mean a and covariance 
matrix C. 

II. System model 

We consider a SIMO channel with R receive antennas. The 
fading in each component SISO channel follows a correlated 
block-fading model 121, with input-output relation for a given 
block 

y™ = Vpdiag(/i„i)a; + n„,, TOe[l:i?] (1) 

where p denotes the signal-to-noise ratio (SNR), x e is 
the transmitted signal vector, y„i G is the received signal 
vector corresponding to the mth receive antenna, and rim ^ 
CN{OJl) is additive noise. Finally, Ka ^ CAf{0,QQ^) is 
the vector of channel coefficients between the transmit antenna 
and the mth receive antenna. HereJl Q g C^^^ and R < Q = 
rank((5) < L. 

Without loss of generality, we assume that the row vectors 
q[ of Q satisfy Qi ^ {i E [1 : L]). The vectors hm and 
n,„ are assumed to be mutually independent and independent 
across m G [1 : i?]. It will turn out to be convenient to write the 
channel-coefficient vector in whitened form as h,,„ = Qsm, 
where s,„ G with Sm ^ CJ\f{0, Iq). Finally, we assume 
that Sjn and n„i change in an independent fashion from block 
to block for all m G [1:R]. 

Setting A (yT . . yT-j^ a (^sJ , . . . , sJ^), = 
{nj , . . . ,1%]^), and X = diag(a;), we can combine the 

^ When Q = L, capacity is known to grow double-Iogaiithmically in SNR 
(7], and, hence, the pre-log is equal to zero. If i? > Q we can always achieve 
the same pre-log as for R = Q by simply using only Q receive antennas. 



individual input-output relations in ([T]i into the overall input- 
output relation 

y = Vpy + with y = (Jh (g) XQ)s. (2) 
III. Lower bound on the pre-log 
The capacity of the channel ^ is defined as 

C{p) = {l/L) sup I {x;y) (3) 

where I{x;y) denotes mutual information 18] p. 251] and the 
supremum is taken over all input distributions f{x) on 
that satisfy the average power constraint E[||a;|p] < L. The 
pre-log is defined as limp^oo(C'(p)/ log(p)). 

The main result of this paper is the following theorem. 
Theorem 1: Suppose that Q satisfies the following 
Property (A): There exists a subset of indices /C C [1 : L] 
with cardinality 

|/C| ^min(r(gi?-l)/(i?-l)l,i) (4) 

such that every Q row vectors of Qic are linearly independent. 
Then, the capacity of the SIMO channel (|2]i can be lower- 
bounded as 

Cip)>[ (l-V^)log(p) + W, if^<L 
- \r{1 - Q/L) log(p) + else. 

Remark 1: The SISO pre-log is 1 - Q/L E). 

Remark 2: For Q — R Theorem [T] reduces to JT] Th. 1]. 

Remark 3: Without loss of generality, we will henceforth 
assume that the set K, in Theorem [T] is given by K, = [1 : |/C|]. 
This can always be achieved by reordering the scalar input- 
output relations in ([T]i. 

Remark 4: The pre-log is (1 — 1/L), provided that L > 
2Q - 1, even if i? = 2 only. 

Remark 5: Property (A) in Theorem [Tjis not restrictive and 
is satisfied for a broad class of correlation matrices. 

Proof: Since we are interested in a capacity lower bound, 
we can evaluate the mutual information in (|3]l for an ap- 
propriate input distribution. Specifically, we take the input 
distribution to have entries Xi (i G [1 : L]) that are independent 
and identically distributed (i.i.d.), zero mean, unit variance, 
and satisfy h{xi) > —oo. This implies that Q Lem. 6.7] 

E[log(|x,|)] > -oo, ie[l:L]. (6) 

For example, we can take Xi ^ CAf{Q, 1). The mutual 
information I{x; y) ~ h{y) — h{y \ x) in ([3]l, evaluated for any 
input distribution satisfying these constraints, is then lower- 
bounded as follows. We first upper-bound h{y \ x) according 
to m Eq. (8)] 

h{y\x)<QR\og{p) + 0{l) (7) 
and then lower-bound h{y) as in Eq. (12)] 

h{y) >iQR + L~a) log(p) + h{Py \xr) + c (8) 



where c is a constant that is independent of p, h{Py \ x-p) is 
independent of p, V ^ [l:a] with a G [1:-^], and 



for sets Xi, . . . , Iji C [1:L] satisfying 
\I^\^QR + L- 



E 

iell:R] 



a. 



(9) 



(10) 



The set V can be interpreted as a set of pilot positions IT]. 
Combining Q and ([8]l, the capacity lower bound in ^ is 
established by choosing 



a 



1, 



if {QR~1)/{R-1) < L 
else, 



QR-{R-l)L, 

provided that we can find sets Ii , . . . , Ij^ such that 
h{Py\x-p) > — oo. The remainder of the paper is devoted 
to identifying such a choice for Ii , . . . , Xji and proving that 
the corresponding differential entropy h{Py \ x-p) is, indeed, 
finite. The main idea is to choose the sets Ii, . . . , I^. such that 
h{Py I x-p) can be related to h{s, x-p) = h{s) + h{x-p,) with 
T) = [q + 1 : L] through a deterministic one-to-one mapping. 
The quantity h{s) + h{x-p) is much easier to deal with than 
h{Py\xv). 

Condition (fTOl i implies that the mapping 

is,xv)^Py^PilB.(E>XQ)s (11) 

is between two vector spaces of the same dimension QR + 
L — a, which is a necessary condition for this mapping to 
be one-to-one. Note that the RHS of ( fTTT i also depends on 
xp, which is, however, taken to be fixed, reflecting the fact 
that the pilot symbols are known to both transmitter and 
receiver Any dependence on x will henceforth implicitly mean 
a dependence on x-p only. We set Ij^ = [1 • ^] and shall choose 
Ii, . . . ,Tr-i C [1:L] as follows: 

(a) If (Qi?-l)/(i?-l) > L, we set Ii = • • • = [1: 
L]. 

(b) If {QR - 1)/{R -1)<L, we let 

2:™ = [l:Q + fc + l], me\l:R-l-l] 
T^ = [l:Q + k + 2], mC,\R-l:R-l] 



(12) 



with k= [{Q-R)/{R-1)\ andl = Q-R-k{R-l) 
Now let 

dPy 



J{s, Xv) 



(13) 



d{s, x-d) 

be the Jacobian of the mapping in ( fTTT i. If this mapping is one- 
to-one on cQ^+i-" almost everywhere (a.e.), we can apply 
the change-of-variables theorem for integrals ||9] Th. 7.26] in 
combination with ifTOl Th. 7.2] and find that 

h{Py I Xp) = h(s,x-v)+2¥.[\og(\Aet{J{s,x-v))\)]. 

The proof is then concluded by establishing that the mapping 
in (fTTT i is one-to-one a.e. and 

E[log(| det(J(s, x-v))\)] > -oo. (14) 

This requires an in-depth analysis of the Jacobian in ( fT3T l. 
which will be carried out in the next section. ■ 



IV. Properties of the Jacobian 

The following lemma provides important insights into the 
structure of the determinant of the Jacobian in (ITjt . 
Lemma 1: The Jacobian in ( fTsT l can be decomposed as 



where 



J{s,x-d) ^ Ji{xd)J2{s)J^{xv) (15) 

Mxr,)^P{lR®X)P'^ (16) 

J2(S) ^P[Jh®Q I fla+l I ••• I fli] (17) 

JaCa;©) = diag(jQfl,diag(a;i5)"^) (18) 



with 



a, ^ (Jfl ® diag(e,)Q)s, i^[l:L\. (19) 
Proof: The lemma follows by noting that 
dy d 



dx; dx 



^ (g) diag(ej)(3)s 

ai, i e [l:L] 



dv 

■jf-=lR(g) XQ. 

OS 

Based on (|6]), we can conclude tha^ 

E[log(|dct(Ji(a;i,))det(J3(a;p))|)] 
= J2 ^'log(l^.l) + E E IE[log(|a;,|)] > -oo. 



jev 



E 

ie[i:i?-i] jei,\v 



To conclude the proof of (fT4l i it therefore remains to show that 
E[log(| det( J2('S))|)] > — oo. Direct computation reveals that 
each vector Pui in ((Tt} with i E Ir \ Ir-i contains only 
one nonzero element, which is given by qj sr. Applying the 
Laplace formula ifTTl p. 7] to J2{s) in ( fTTI i therefore yields the 
decomposition 

|dct(J2(s))| = |dct(J4(s))| J] \qJsR\ (20) 

i6lfi\Ifi_i 



with 



and 



Ji{s)^Pi[{lR®Q) I a„+i 



(21) 



Pi = diag((/L)ii, . . . , 

The expectation of the logarithm of the second term on the 
RHS in ( I20I 1 is finite because sr ~ CAf{0, Iq). It remains 
to show that E[log(| dct(J4(s))|)] > —00. This is the most 
technical part in the proof of Theorem [T] and can be accom- 
plished by applying methods from algebraic geometry, namely 
Theorem 12] in Appendix lAl which is a consequence of Hiron- 
aka's celebrated theorem on resolution of singularities H, ||5|. 
A direct proof would require showing that the expected log 
of the determinant of the (high-dimensional) matrix J4(s) is 
finite, which seems exceedingly difficult. Hironaka's theorem 

^^We assume that xj 7^ for all j £ V. 



drastically simplifies the proof as it teUs us that dct(J4(-)) ^ 
implies that E[log(| dct(J4(s))|)] > — oo. We start by noting 
that det(J4(s)) is a homogeneous polynomial in si, . . . , sqr 
of degree D = \Ir-i \ — a, i.e.. 



If E 



det(J4(As)) = det(J4(s)), VA G 



(22) 



which allows us to apply the following proposition (for AI = 
QR and D defined above): 

Proposition 1: Let 5 be a homogeneous polynomial in 
si,...,SM of degree D e No with s ^ CN{0,Im)- Then 
g ^ implies that 



E[log(|(7(s)|)] > -00. 



(23) 



Proof: Writing s as s = and using the fact 

that g is a homogeneous polynomial of degree D, we can 
upper-bound the absolute value of the expectation in (|23] | by 

|E[log(|g(^)|)]| <i?E[|log(||^||)|]+E[|log(|g(s/||s||)|)|]. 



Then ||7] Lem. 6.7] together with s ^ CA/^(0, Im) implies that 
A < 00. Introducing polar coordinates lfT2l p. 55] r e and 



= (^^l,. 



J2M 



2A/-2 



X [0, 27r] 



for the complex vector s G C^^, we can further upper-bound 
B according to 

B< ( cxp(-r2)r2^^-idr X j | log(|/(0)|)| d6> 

JO J A 

where / is obtained from g by changing to polar coordinates 
r G € A. Note that / is a real analytic function ||6] 

Def. 1.6.1], 5^0 implies that / ^ 0, and we are integrating 
|log(|/(-)|)| over a compact set A. We can therefore apply 
Theorem |2] in Appendix lAl to conclude that B < 00. ■ 

It remains to show that dct(J4(-)) ^ for our specific 
choice of sets Ii, . . . which will be proved in the 

following lemma: 

Lemma 2: Property (A) in Theorem [T] implies that 
det(J4(-)) ^ 0. 

Proof: See Appendix IbI ■ 

In summary, we proved that (fl4l holds provided that 
Property (A) in Theorem [T] is satisfied. It turns out that 
det(J4(-)) ^ also implies that the mapping in (fTTI ) is one- 



to-one a.e. on 



The proof, which is along the lines 
of the proof of |[T] Lem. 2], is omitted due to space limitations. 
This completes the proof of Theorem [T] 

Appendix 
A. Resolution of singularities 

In this appendix, we show how Hironaka's theorem 
on resolution of singularities can be used to prove that 
/acR" Uos(I/(^)I)M^ < °° provided that / ^ is a real 
analytic function and A is a compact set. 

We start by defining notation that will be used in this 
appendix. Let C^{x,e) = {y e R^^ | \xi - yi\ < e, V i e 
[1:A/]} denote the open cube with side length 2e centered at 



X. For X e C^^ and m e N^^ let x"^ = ^ 
is a subset of the image of a map / then denotes the 

inverse image of S. 

The following lemma is an immediate consequence of a 
modified version of Hironaka's theorem ||3] Th. 2.3]. This 
modified version originally appeared in [[T3l . The main point 
of this lemma is that it allows us to rewrite any real analytic 
function ||6] Def. 1.6.1] / ^ locally as a product of a 
monomial and a nonvanishing real analytic function. 

Lemma 3: Let / be a real analytic function from a neigh- 
borhood 17 C of a; e to R. Suppose that f{x) = 0. 
Then, there exists a triple (S, 5, c^), where 

(a) S C is an open set in M.^ with x E Y,, 

(b) S is a ZJ-dimensional real analytic manifold ^ Def. 2.10], 

(c) : S — > E is a real analytic map 

that satisfies the following conditions: 

(1) The map (p is proper, i.e., the inverse image of any compact 
set is compact. 

(2) The map is a real analytic isomorphism between S \ 
(/o </,)-! (0) andl]\.ri(0). 

(3) For each point P E S, there exists a coordinate chart 
{Ep,(pp} such that P e Ep C E, ipp : C^{0,ep) ^ 

is a real analytic isomorphism for some ep > with 
(^p(O) = P, 

\f{(Poipp{u))\ = /ip(m)m"^^, Vm G C^(0,ep) 

where hp is a nonvanishing real analytic function on 
C^(0, ep) and mp G Nq, and the determinant of the 
Jacobian of the mapping (po ipp satisfies 

det (^(^^|£M) = v« e C-{0,ep) 

where gp is a nonvanishing real analytic function on 
C^(0,ep) and np e N^. 

Proof: The main idea is to apply |3l Th. 2.3] to the 
function f{z) = f{z + x), Vz E il — x. We omit the details 
due to space limitations. ■ 

We are now in a position to state the theorem that is needed 
to prove i? < 00 in the proof of Proposition [T| 

Theorem 2: Let / ^ be a real analytic function on an 
open set C M^. Then 



|log(|/(a;)|)|da;<oo 



(24) 



for all compact sets A C fi. 

Proof: Let a; G A. If f{x) = then Lemma [3] implies 
that there exists a triple (Sg;, S3;, 0a;) where S^; C is an 
open set containing x, is a real analytic manifold, and 
(j)x :Ex Sa; is a proper real analytic map. Moreover, for 
each P E Ex there exists a coordinate chart {E^.p, (px,p} 
such that Ex,p = 'fia:,p{C^{0,ex,p)), Vx,p(0) = P, and 



det 



\f o (j)xO (px^p{u)\ = hx,p{u)u^'^ '' 



du 



gx,p{u)uJ 



(25) 



for all u e C^(0,ea;.p), where g^ p and h^^p are non- 
vanishing real analytic functions on {O^Cx^p). We can 
choose e^ p sufficiently small so that g^ p and hx,p are 
bounded on C^(0, ex,p)- If .f{x) ^ the existence of a triple 
(S^, Sj., with the properties specified above is guaranteed 
by taking S^; = Sa, sufficiently small such that / does not 
vanish on and by setting (^x to be the identity map. 
Now for each a; e A, we choose an open neighborhood 
and a compact neighborhood A^, such that a; e C 
Aa; C Sa;- Sincc A is a compact set, there exists a finite set 
of vectors {x\, . . . , a;jv} in A such that 



A c U 

i6[l:Ar] 



c 



u 

ie[l:Af] 



„ A — 



For each i e [1 : N\ set A, = A^^,, = Y^x,. 
and 0i = 0a;i ■ Since the mapping (^i : S.; — > is proper, each 
set 0i~^(Ai) C Si is a compact set. Therefore, there exists a 
finite set of points {Pi, ... , Pa/.} in such that 



je[l:A/,] 



(26) 



with ^i^j = .^aj. p.. Since (|26] l holds for all i e [1 : A^], we 
can upper-bound the integral in (l24l i as follows: 

log(|/(a;)|)|da;< ^ / | log(|/(a;)|)| da; 

ie[l:Af]"^'^' 



E E / |M"'-l0g(|M-'-|)|dM + C2 

< oo 

where ci,C2 > are positive real numbers, gij,hij are 
bounded nonvanishing real analytic functions on C^(0, ei.j), 
rriij, riij are vectors of nonnegative integers, and we changed 
variables according to (|25t . ■ 

B. Proof of Lemma |2] 

We present a proof for a ^ 1,1 ^ and skip the (simpler) 
cases a > 1 and a = 1, / = 0. 

Suppose that a = 1 and / ^ 0. We can write J4{s) in (fTTI i 
as J4(s) = [diag(Qi,,...,Qi^_i,Qx^_J A] with ^ 
[Aj'' . . . A]^] and Am defined as 



/ 

T 
92 «r: 



V 



/ 
T 
92 



V 



T 



0\ 




0/ 



me [1:P-/-1] 



T 

lQ+k+2' 



me[R-l:R] 



Property (A) in Theorem [T] implies that for arbitrary subsets 

/Cm C [2:Q + fc + l], me[l:R-l-l] 
/C™ C [2:Q + fc + 2], me[R-l:R] 

with |/Cm| = Q— 1, we can find vectors s™ G (m e [1 :i?]) 
such that 

(a) qjs„i = for all vectors qj with j E /Cm; 

(b) qij Sm 7^ for all vectors with j <E IC\ Km- 

This implies that for each choice of such sets /C,„ (m G [1 : i?]), 
there exists a set of vectors s„i G (m G [1 : i?]) such 
that the number of nonzero elements rim in each matrix Am 
satisfies 



(27) 



'^m = fc + l, mG[l:i?-/-l] 
= + 2, m e[R- l:R]. 

Moreover, we have 

n,„ = Q + fc + 1 

me[l:_R] 

which implies that we can choose the subsets JCm and the 
vectors Sm (ni G [1 : R]) such that each column of A contains 
precisely one nonzero element. Applying the Laplace formula 
ifTTl p. 7] iteratively, we therefore get 

I det(J4(s))| = c Jl I det(Q;c„u{i})l > 

me[l:-R] 

where c is a positive constant and we used Property (A) in 
Theorem [U in the last step. ■ 
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